Spin-lattice relaxation in the mixed state of the high- Tc cuprates: 
electronic spin-flip scattering versus spin-fluctuations 
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Recent experimental and theoretical studies have established that the spin-lattice relaxation rate, 
1/Ti, measured in nuclear magnetic resonance (NMR) experiments is a site-sensitive probe for the 
electronic spectrum in the mixed state of the high- Tc cuprates. While some groups suggested that 
l/Ti is solely determined by spin-flip scattering of BCS-like electrons, other groups stressed the 
importance of antiferromagnetic spin-fluctuations. We show that these two relaxation mechanisms 
give rise to a qualitatively different temperature and frequency dependence of ^^O l/Ti. A compari- 
son of our results with the experimental ^'^O l/Ti data provides support for a relaxation mechanism 
dominated by antiferromagnetic spin-fluctuations. 

PACS numbers: 74.25.Nf, 74.60.Ec, 74.25.Ha 

Recently the temperature and frequency dependence of 
the ^^Cu and ^^O spin-lattice relaxation rate, 1/Ti, in the 
mixed state of YBa2Cu306+2; (YBCO) has been the focus 
of intense experimental efforts . It was first shown by 
Curro et al. (CMHS) and subsequently by Halperin 
et al. that ^^O 1/Ti in the mixed state of YBa2Cu307 
increases with increasing resonance frequency, i.e., with 
decreasing distance from the vortex core. Preliminary 
results of Milling et al. ^ show a qualitatively similar 
behavior for the ^^Cu relaxation rate in YBCO. CMHS 
were the first to point out that this frequency dependence 
of 1/Ti arises from the presence of a supercurrent, which 
gives rise to a site-specific density of states, and conse- 
quently a site-specific 1 /Ti . This interpretation receives 
strong support from the experimental observation that 
the frequency dependence of 1/Ti disappears above the 
melting temperature of the vortex lattice. 

Several theoretical models [^-Q have been suggested 
to explain the experimentally observed frequency depen- 
dence of 1/Ti. While it was argued in Refs. that 
1/Ti for ^•^Cu arises solely from electronic spin- flip scat- 
tering of BCS-type electrons, Morr and Wortis (MW) ||] 
proposed that antiferromagnetic spin-fluctuations pro- 
vide the dominant contribution to the ^'^Cu spin-lattice 
relaxation rate. MW also conjectured that the combined 
temperature and frequency dependence of 1/Ti can dis- 
tinguish between relaxation due to electronic spin-flip 
(ESF) scattering and relaxation due to antiferromagnetic 
spin- fluctuations (ASF) . 

In this communication we show that due to the spe- 
cific form of the "'^^O hyperfine coupling, the tempera- 
ture and frequency dependence of the ^''O spin-lattice 
relaxation rate for ESF and ASF relaxation is qualita- 
tively different. In particular, we demonstrate that for 
the ESF mechanism, l/Ti increases monotonically with 
temperature and possesses a large distribution of values, 
for a given resonance frequency, close to the vortex core. 
In contrast, for the ASF relaxation mechanism, we ob- 
tain a well-defined relation between l/Ti and resonance 
frequency, and a non-monotonic dependence of 1/Ti on 
temperature. Though a comparison of our theoretical re- 



sults with the available experimental data suggests that 
the spin-lattice relaxation rate in the mixed state is dom- 
inated by the ASF mechanism, further experiments are 
required to test the predictions presented in this paper. 
To the extent that future experiments confirm our con- 
clusion, it immediately implies that the relaxation rate 
of ^■^Cu is as well dominated by antiferromagnetic spin 
fluctuations due to its even more favorable hyperfine cou- 
pling. 

We begin by shortly reviewing the theoretical models 
for an ESF and ASF based spin-lattice relaxation rate; a 
detailed account is given in Refs. For ESF relax- 

ation 1^ , we consider the model of Ref. Q , whose results 
are in good qualitative with those presented in Ref. Q . 

The general expression for the ^^O spin-lattice relax- 
ation rate, 1/Ti, with an applied field parallel to the 
c— axis is given by 
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where 
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^^0(2) (q) = cos(gj2)2 Fqo) (q) = cos(g,/2)2 , 

(2) 

C is the ^^O transferred hyperfine coupling constant, and 
0(2) [0(3)] refers to the oxygen nucleus which is located 
between two Cu nuclei along the x(y)-axis. 

Within the ESF approach the spin-susceptibility in the 
superconducting state, neglecting final state effects, is 
given by % = E, where 



S(q, iu,,) - -T ^ {G(k, ir!„OG(k + q, iVL,, 
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and G and F are the normal and anomalous Green's func- 
tions 
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In the limit that the supercurrent momentum, p^, varies 
on a length scale larger than 1/kp, the dispersion of the 
fermionic quasi-particles, £^k, up to linear order in ps is 
in semiclassical approximation given by M 
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where et is the electronic normal state dispersion, Ak = 
Ao (cos(fcj;) — cos(A:y))/2 is the d-wave gap, and Vi?(k) = 
9ek/9k. Following MW, we neglect the electronic 
Zeeman-splitting which for typical applied fields is much 
smaller than the second term in Eq.(^, the so-called 
Doppler-shift, as well as any effects of the supercurrent 
on Ak, Uk and Vk 0. 

Spin relaxation due to the ASF mechanism is described 
by the spin-fermion model ||ll| where the spin susceptibil- 
ity, Xj is renormalized by the interaction with low-energy 
fermionic degrees of freedom and given by 
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Here 
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is the bare propagator, is the bare magnetic correlation 
length, a is a temperature independent constant, Q = 
(tt, tt) is the position of the magnetic peak in momentum 
space [0 and 11 is the bosonic self-energy given by the 
irreducible particle-hole bubble. A supercurrent affects 
X only through H, which is calculated using second order 
perturbation theory in the spin-fermion coupling, g, and 
one obtains 
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with E in Eq.(H). In the limit lu — > 0, Re 11 is only 
weakly momentum dependent and can thus be included 
in a renormalized but momentum independent correla- 
tion length, ^, in Eq.(0). 

For both the ESF and ASF mechanism, we are thus 
left with the calculation of Im S which for is dom- 
inated by particle-hole excitations in the vicinity of the 
superconducting nodes. Expanding i?k up to linear or- 
der in momentum around the nodes, we can perform the 
momentum and frequency integrations in Eq.(^) analyt- 
ically. Combining the resulting expression for Im S with 
Eqs.(^, (^) and (|^) we obtain the spin-lattice relaxation 
rate presented in Eq.( [Al| ) of the appendix. For our sub- 
sequent discussion of the ESF and ASF mechanism it is 
necessary to briefly review the form of 1/Ti in two limits. 
A low temperatur es w here \Dm/T\ 3> 1 for m = 1..4, one 
obtains from Eq.(|All) 
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where J^{qn,m), Dm are defined in the appendix, and the 
sum runs only over those nodes with D^n-, -D„ < 0. Since 
the supercurrent in general breaks the lattice symmetry, 
the 0(2) and 0(3) relaxation rates are expected to be 
different, as we shall discuss below. In contrast, in the 
high-temperature limit {\Dm/T\ <^ 1 for m — 1..4) one 
has 



T 



1 _ttA f 
T\f " 3N \VFVA 



1 4 

E 

n.m—l 



^(qn 



(10) 



and we recover the temperature dependence of the relax- 
ation rate in the absence of a supercurrent. 

A site specific relation between the resonance fre- 
quency A:^(r) — 7„?iiJ;oc(r) where 7„ is the ^''O nuclear 
gyromagnetic ratio and Hioc{y) is the local magnetic field, 
and 1/Ti is obtained via the local supercurrent momen- 
tum Ps(r) V X H;oc(r); a relation which is valid even in 
the presence of non-local as well as non-linear corrections 
]l3[ . Following MW, we restrict our discussion of the re- 
laxation rates to nuclei which are further than R > 2^ab 
from the center of the vortex core (where Cab is the super- 
conducting in-plane coherence length) Q. We consider 
a hexagonal vortex lattice [ p^p6| and use the parameter 
set VF « 0.4 eV, va « 20 meV, and ^ « 2 0, which was 
shown to describe fermionic and magnetic excitations in 
YBaaCusOy |Jl|. 

In Fig. |1| we plot the ^^O spectrum H which describes 
the distribution of local magnetic fields (solid line) and 
1^0(2,3) 1/TiT (open squares) for (a) the ASF mecha- 
nism, and (b) the ESF mechanism as a function of Hioc, 
i.e., resonance frequency. Nuclei at the highest frequen- 
cies are located at a distance 2^ab from the center of the 
vortex (nuclei closer to the vortex core have been omit- 
ted, see above discussion), nuclei at the lowest frequencies 
are in the center of a vortex triangle, those at the max- 
imum of the spectrum are at the midpoint between two 
vortices. 

For both relaxation mechanisms, 1/TiT increases with 
increasing resonance frequency, i.e., decreasing distance 
from the vortex core, similar to the results obtained in 
Refs. for ^^Cu. For T = 1 K, the relaxation rate 

for practically all nuclei is given by l/TiT ~ p^, Eq.(||), 
and since |ps| increases with decreasing distance from 
the vortex core, l/TiT increases consequently. However, 
l/TiT for the ESF mechanism exhibits a large distri- 
bution of values for a given Hioc, in contrast to l/TiT 
for ASF relaxation. Nuclei with the same resonance fre- 
quency can in general possess different relaxation rates, 
since 1/TiT depends not only on the magnitude of p^, 
but also on the angle between the direction of the local 
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FIG. 1. ^''O spectrum (solid line) together with ^^0(2,3) 



(l/TiT) as a function of local magnetic field, Hioc, 
T = IK. (a) ASF mechanism, and (b) ESF mechanism. 



for 



supercurrent and the crystal axes, i.e., the Fermi veloc- 
ity vp at the nodal points, Eq.(|9|). It is exactly this 
angular dependence which for ESF relaxation leads to 
the large distribution of values for l/TiT. The absence 
of such a large distribution for the ASF mechanism is 
due to the antiferromagnetic enhancement factor (i.e., 



the Stoncr enhancement) in the denominator of Eq.(A3). 
The Stoner enhancement favors particle-hole excitations 
which connect nodes at opposite sites of the Fermi sur- 
face (FS), i.e., with large momentum transfer, qn,m = 



(see Fig. ^). Due to this enhancement, the second term 
on the r.h.s. of Eq.(|^) yields the dominant contribution to 
1/TiT for ASF relaxation in the low-temperature limit, 
and one obtains 
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Thus the leading temperature contribution to l/TiT is 
independent of the angle between Ps and the underly- 
ing lattice, and the narrow distribution of 1/TiT seen at 
higher frequencies for ASF relaxation arises solely from 
subleading contributions including particle-hole excita- 
tions with small momentum transfer, e.g., with qn m = k^. 
(see Fig. |a). 
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FIG. 2. (a) Fermi Surface and particle-hole excitations 
(dashed arrow) with q^.m = k;,a (nodes of the superconduct- 
ing gap are shown as dotted lines), (b) Angle between ps and 
the crystal x-axis. 

In Fig. H we present the ^^0(2,3) relaxation rate for 
three different temperatures (while the 30 K and 60 K 
curves are horizontally offset for clarity, the left end 
points of all curves correspond to the same local field 
at a distance 2£^ab from the vortex core). The tempera- 
ture dependence of 1/TiT in the ASF approach (Fig. Ha) 
is similar to that of ^^Cu shown in Fig. 2 of Ref. [pj. 
With increasing temperature, 1/TiT at high resonance 
frequencies decreases, in agreement with Eq.(pT|), while 
at low frequencies 1/TiT ~ T^, Eq.([Tc|), and the relax- 
ation rate increases. As a result, 1/TiT possesses a non- 
monotonic dependence on the local field, as shown by the 
minimum in the relaxation rate for intermediate frequen- 
cies, which was discussed in detail in Ref. 

In contrast, 1/TiT for the ESF mechanism increases 
with increasing temperature at all resonance frequencies. 
At low frequencies, i.e., in the high temperature limit, 
this follows directly from Eq.(lO). At higher frequen- 
cies, i.e., in the low temperature limit, an analysis of 
Eq.(H) shows that due to the ^^O form factors, the first 
term on the r.h.s. of Eq.(^ dominates 1/TiT which there- 
fore also increases with increasing temperature. Conse- 
quently, 1/TiT does noi exhibit a minimum as a function 
of the local field as was the case for ASF relaxation. Note 
also that the distribution of values for l/TiT decreases as 
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FIG. 3. ^■'0(2,3) (l/TiT) as a function of local field for 
three different temperatures, (a) ASF mechanism, (b) ESF 
mechanism. The 30 K and 60 K curves are horizontally offset 
for clarity. 

the temperature increases, since in the high temperature 
limit the relaxation rate becomes independent of reso- 
nance frequency, i.e., 1/TiT T^. The spin relaxation 
rate thus possesses qualitatively different frequency and 
temperature dependencies for the ASF and ESF mech- 
anism, the origin of which lies in the antiferromagnetic 
enhancement factor in the denominator of Eq. ( |A3| ) which 
is only present for ASF relaxation. 

We next compare our theoretical results with the ex- 
perimental data on i^O(2,3) 1/TiT by CMHS ||, which 
we present in Fig. |^. The experimentally measured 
1/TiT increases with increasing resonance frequency, in 
agreement with our theoretical results in Fig. ^, decreases 
between T=5K and T=10K with a larger reduction at 
higher frequencies, and increases between T=10K and 
T=40K. No evidence for a distribution of relaxation rates 
at high frequencies was found so far, however, a fur- 
ther quantitative analysis is still required to determine 
an upper bound for the spread in l/Ti |p^ . Though 
both experimental features agree qualitatively with our 
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FIG. 4. Experimental data for the ^^O spectrum (solid 
line) together with (l/TiT) as a function of resonance fre- 
quency for three different temperatures. Data are taken from 



theoretical predictions for I/TiT in the ASF approach 
(see Fig. spin-diffusion, impurity effects and vortex 
vibrations |l9|] can potentially contribute to the spin- 
lattice relaxation and thus complicate the analysis of 
the experimental data within the ASF or ESF approach. 
While it was argued that spin-diffusion is strongly 
suppressed by the inhomogeneity of the magnetic field 
in a vortex lattice, CMHS concluded from the different 
temperature dependence of the apical and planar oxy- 
gen relaxation rates, that for T > 25 vortex vibrations 
are irrelevant for the relaxation of ^^O; thus the relax- 
ation arises solely from the electronic/magnetic mecha- 
nisms discussed above. For T < 25, the apical and pla- 
nar oxygen relaxation rates exhibit the same linear tem- 
perature dependence, however, the planar I/Ti is still 
considerably larger than the apical l/Ti. Though this 
result does not exclude the possibility that vortex vibra- 
tions contribute to the relaxation process below 25 K, 
it suggests that the relaxation mechanism is still dom- 
inated by electronic/magnetic excitations. Clearly, fur- 
ther measurements of 1/TiT are required to study the 
various relaxation mechanisms pO[ |. However, to the ex- 
tent that the effects of vortex vibrations and impurities 
are negligible at low temperatures, the experimental data 
in Fig. 1^ support a predominant ASF mechanism of the 
spin relaxation in the mixed state. 

Since the supercurrent in general breaks the symme- 
try of the underlying lattice, one expects different relax- 
ation rates for 0(2) and 0(3) for a given direction of 
Ps. An experiment which is able to reveal this differ- 
ence in the relaxation rates is a nuclear quadrupole res- 
onance (NQR) measurement in which the direction of a 
uniform supercurrent is varied. In Fig. |^a we plot l/TiT 
for the ASF mechanism as a function of the angle, (j>, 
between the supercurrent momentum, p^, and the crys- 
tal axes (see Fig. ||d). It follows from Eq.(^ that in the 
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FIG. 5. "0(2,3)l/rir for ps = 0.02 and three different 
temperatures as a function of the angle (j> between and the 
crystal x— axis (see Fig. Ijb). (a) ASF mechanism, (b) ESF 
mechanism. 



limit T ^ 0, 0(2) and 0(3) 1/TiT possess the same an- 
gular dependence. This identical behavior again results 
from the Stoner enhancement, since for the dominating 
particle- hole excitations with large momentum transfer, 
q_n,m = k;, (sec Fig. |), the form factors F{qn,m) for 
0(2) and 0(3) are the same. With increasing temper- 
ature, the relaxation rates for 0(2) and 0(3) begin to 
deviate, since excitations with smaller momentum trans- 
fer become increasingly more important, however, the 
maxima and minima of 0(2) and 0(3) \/TiT still occur 
at the same angle 4>. 

The relaxation rates of 0(2,3) for the ESF mechanism, 
(Fig. ^), are also identical in the low-temperature limit, 
as follows directly from Eq.(|^). However, with increasing 
temperature, the angular dependence of 1/TiT for 0(2) 
and 0(3) is out of phase: the minima of 0(2) \/TiT co- 
incide with the maxima of 0(3) 1/TiT. To demonstrate 
that this "phase-shift" originates from the different form 
factors of 0(2,3), Eq.(^, we consider the case (j) = n. 
One then has from Eq.(S) for 0(2) 
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whereas for 0(3) one obtains 



1 + cos'{k^/2)]{Dl + Dl + 2D1D4) + 0{T^) . (13) 
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For T = we indeed find that the 0(2,3) relaxation rates 
are identical, while for T > (and (j) = it) 0(2) l/TiT is 
larger than 0(3)l/TiT, in agreement with our numerical 
results in Fig. ||. 

Finally, while the overall scale of l/TiT within the ESF 
and ASF approach depends on a variety of parameters, 
e.g., the hyperfine coupling, C, the magnetic correlations 
length, ^, the detailed form of the Fermi surface, etc., 
we find that its qualitative features are robust against 
changes in these parameters. Thus our predictions for 
the qualitative differences between the ESF and ASF re- 
laxation mechanisms are valid for all cuprate materials, 
and do not depend on any fine-tuning of parameters. 

In summary, we present a theoretical scenario for the 
spin-lattice relaxation rate of ^^O in the mixed state of 
the high-Tc cuprates. We consider 1/Ti for two relax- 
ation mechanisms: one based on electronic spin-flip scat- 
tering of BCS-type electrons, and one due to antiferro- 
magnetic spin fluctuations. We show that the predicted 
temperature and frequency dependence of 1/Ti differs 
qualitatively for these two mechanisms. A comparison 
of our theoretical results with the available experimen- 
tal data from Ref. ||^ suggests that the ASF mechanism 
dominates the relaxation rate. To the extent that this 
conclusion is confirmed by future more detailed experi- 
ments, it implies that the relaxation rate of ^^Cu in the 
mixed state is also dominated by the ASF mechanism, in 
contrast to the scenarios pursued in Refs. Finally, 
we propose an NQR experiment in which the direction 
of a uniform supercurrent with respect to the crystal lat- 
tice is varied and predict a unique angular dependence of 
1/Ti for 0(2) and 0(3). 
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National Laboratory. 



APPENDIX A: 

The general form for l/TiT in the presence of a super- 
current is given by 
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where Dm, = ■ p^, is the Fermi velocity at 

the node in the m'th quadrant of the BriUouin zone, 
va — |9Ak/9k| at the nodes, and qn,m is the momentum 
connecting nodes n and m. For ESF scattering we find 



while for the ASF mechanism on has 

A^{ag)^kB{n^jnlef/{2n) 

-Pb(2,3)(qn,m) 
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In both cases, Gi_2 is given by 
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where np{e) = [exp(e/T) + 1] ^ is the Fermi function, 
e = max{Dn, D„i}, and 
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